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ABSTRACT 


We present a multi-wavelength observational study of a low-mass star-forming 
region, L1251-C, with observational results at wavelengths from the near-infrared 
to the millimeter. Spitzer Space Telescope observations confirmed that IRAS 
22343-1-7501 is a small group of protostellar objects. The extended emission to 
east-west direction with its intensity peak at the center of L1251A has been de¬ 
tected at 350 and 850 fj,m with the CSO and JCMT telescopes, tracing dense en¬ 
velope materials around L1251A. The single-dish data from the KVN and TRAO 
telescopes show inconsistencies between the intensity peaks of several molecular 
line emission and that of the continuum emission, suggesting complex distribu¬ 
tions of molecular abundances around L1251A. The SMA interferometer data, 
however, show intensity peaks of CO 2-1 and 2-1 located at the position of 
IRS 1, which is both the brightest source in IRAC image and the weakest source 
in the 1.3 mm dust continuum map. IRS 1 is the strongest candidate for the 
driving source of the newly detected compact CO 2-1 outflow. Over the whole 
region (14'xl4') of L1251-C, 3 Class I and 16 Class II sources have been detected, 
including three YSOs in L1251A. A comparison with the average projected dis¬ 
tance among 19 YSOs in L1251-C and that among 3 YSOs in L1251A suggests 
L1251-C is an example of low-mass cluster formation, where protostellar objects 
are forming in a small group. 

Subject headings: ISM; individual objects (L1251-C) - ISM; jets and outflows - stars; 
formation - stars; protostars 
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1. INTRODUCTION 

Many studies of low-mass star formation have been done in details. Star formation 


mostlv occurs in groups and clusters w 

lere 

the number density of young stellar objects 

(YSOs) at least ~3 pc ^ ( 

Lada & Lada 

2003 

). The Spitzer legacy proiect “From Molecular 

Cores to Planets Forming Disks” (c2d; 

Evans et al. 

2003, 

2009 

) found that most sites of 


low-mass star formation within 400 pc have signihcant populations of protostars forming in 
small groups of multiple objects (e.g., <10 members) located in close proximity. 

L1251 is a s tar-forming region located in Cepheus at a distance of 300 ± 50 pc 


flKun et ah 


19931) from the Sun. There are hve dense cores A to E 


cm ^) in LI251, which have been assigned A to E with Right Ascension fISato et al 


199411 


Two outflows, an extended CO 1-0 outflow and a well- collimated compact CO 1-0 outflow 


were previously known in cores C and E, respectiveh 

/ (Sato & Fukui 1989). These two 

outhows are named LI251A and LI25! 

C^®0 core among hve cores in L1251 ( 

.B bv Hodann 

1 

U9941). Core C is the second densest 

Sato et al 

199 

4) where the strongest extended CO 

1-0 outhow lL125lAj is located (Sato & Fukui 

1989 

). The dense region around IRAS 


22376-1-7455, which is the driving source of the we 

1-collimated compact outhow (L1251B) in 

core E has been studied in detail ( 

Lee et al. 

20 

06, 

2007 

). Recently, another well-collimated 

outhow in core A has been revealed ( 

Lee et al. 

201o|. In this paper, we use “L1251-C” to 


present the extended region observed in C^®0 1-0 and it is the same region with “L1251C” 
in the Spitzer data archive. On the other hand, “L1251A” is used to refer the densest region 
around IRAS 22343-1-7501, which is believed to be the driving source of the extended CO 
2-1 outflow. Our Spitzer observations revealed that L1251A is a multi-protostellar system 
consisting of three sources. 


Earlier studies already showed that 


protostellar objects. 


RAS 22343-1-7501 is composed of multiple 


Rosvick fc Davidgel fjl995l ) found that there are hve near-infrared 
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(NI R) sources labe 


and 


Reipurth et al 


ed A- 


around L1251A. Meehan et al. ( 1998 ). Beltran et al. ( 2001 ). 


(120041) presented also mnltiple radio continnum sonrces, detected 


with Very Long Array (VLA), in L 


(iBeltran et ah 


200 


Rosvick fc Davidgel fjl995l ). respecti vely. 


Reipurth et ah 


251A . Two VLA radio sources, VLA 6 and VLA 7 


photometry for sources A and B of 
as classical T Tauri stars. 


2004) are assoc i ated w ith the NIR sources, D and A of 


Knn et al. 


(120091) pr esented optical spectra and 


Rosvick fc Davided (119951) and classified both sources 


N H.'^ cores were also de tected around L1251A (lAndada et ahl 119971: iToth &: Walmslev 


19961) . 


Nikolic et al. 


( 2003 ) presented t 


le distributions o 


compared with the NH 3 distribution of 


Toth fc Walmslev 


several molecular transitions 


(119961) and showed outflow-like 


structures, which ar e possibly driven 


were also detected (Toth fc Walmslev 


Claussen et al. 


jy VLA 6 and VLA 7, in HCO^ 1-0. H 9 O mase rs 


19941: 


Wilking et al. 


1994 : 


Xiang fc Turner 


1995 


19961) in this region. 


Schwartz et al. 


( 1988 ) discovered a bipolar CO 1-0 outflow in L1251A. This CO 


1-0 bipolar on 


flow shows asymmetric blue and red lobes in morphology and intensity. 


Sato fc Fukuil (119891) also detected an extended C O 1-0 outflow in L125 1A with symmetric 


structure. Interestingly, the CO 1-0 outflow from 
southwest (N E-SW) axis, wh i ch is perpendicular to the northwes 


Sato fc Fukuil (119891) has a northeast 


axis found by 


Schwartz et al 


(119881) . An optical jet was found by 


-southeast 


Balazs et al 


( NW-S E) 

( 1992 ) near 


L1251A, supporting L1251A as the origin of the extended CO 1-0 outflow. 


In this paper, we analyse the properties of young stellar objects (YSOs) based on the 
c 2 d survey data and newly obtained molecular line data to study the environment of star 
formation in L1251-C, especially around L1251A. In this study we use the {a,6) coordinates 
of IRAS 22343+7501 as (22^35"*23k4, +75°17'07") in the J2000.0 epoch. 


Our observations of L1251-C are summarised in section 2. Results and analysis of the 
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data are described in section 3 and section 4, respectively. In section 5, we discuss star 
formation in small groups and outflows in L1251A. Finally, a summary of this paper is 
given in section 6. 


2. OBSERVATIONS AND DATA 
2.1. Spitzer Space Telescope Data 


In the Spitzer legacy program c2d, L12 51-C was observed at 3.6, 4.5, 5.8, and 8.0 
/rm with the Infrared Array Camera (IRAC; 


Fazio et ah 


20041) on 2004 October 18 and 


also ob served at 24 and 70 pm. with the Multi-band Imaging Photometer for Spitzer 


(MIPS; 


Rieke et al 


2004) on 2004 September 24. The High-Dynamic Range (HDR) mode 


observations were obtained to get more accurate intensity levels of saturated sources. In 
the HDR mode, a long-integration frame of 10.4 second and a short-integration frame of 
0.4 second were taken at each position. The descripti on of the data proc essing details can 


be found in the c2d hnal delivery documentation (see 


Evans et al. 


2007 )h In this paper. 


we use the c2d photometric data catalog of YSOs including Two Micron All Sky Survey 
(2MASS) photometric data. For the saturated sources, however, we use the photometry 
obtained with the HDR mode data. 


^see http://irsa. ipac.caltech.edu/data/SPITZER/docs/spitzermission/observingprograms/legacy/ 
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2.2. Radio Single-dish Observations and Data 

2.2.1. Korean VLB I Network Observations 

Molecular line observations of L1251-C were carried out with the Korean VLBI Network 
(KVN) 21 m radio telescope at the Yonsei site, Seoul from 2012 October to 2013 March and 
at the Ulsan site, Ulsan in 2013 April. The central position of mapping observations was 
the position of IRAS 22343+7501. We observed the HCO+ 1-0, 1-0, N 2 H+ 1-0, 

and HCN 1-0 lines to map the cloud structure. We also observed the H 2 O 6i6-523 maser 
and the H 2 CO 2i2-lii transition simultaneously. The KVN radio telescopes were designed 
for simultaneous observations in the 22, 43, 86, and 129 GHz bands (Han et ah 2008). The 
pointing was checked every 2 hours using R Gas, a strong SiO maser source near the target 
source. 

The system temperatures during the observations ranged from 60 K to 250 K (at 22 
GHz), from 100 K to 250 K (at 86 GHz), and from 150 K to 400 K (at 129 GHz), depending 
on weather conditions and elevations. We used digital spectrometer in its 32 MHz mode, 
yielding velocity resolutions of 0.21 km s“^ (4096 channels at 22 GHz), 0.21 km s“^ (4096 
channels at 86 GHz), and 0.27 km s“^ (4096 channels at 129 GHz). The 1 a rms noise levels 
are presented in Table 1. The data were calibrated using the chopper wheel method to get 
an antenna temperature T^. 

The beam Full Width at Half Maximum (FWHM) is ~30 " at 86 GHz. The maps were 
sampled every a half of beam FWHM (15"). We used the GLAS^ reduction package from 
Institue de Radioastronomi Millimetrique (IRAM) for the data reduction. A summary of 
observations is presented in Table 1. 


^see http://www.iram.fr/IRAMFR/GILDAS/doc/html/class-html/class.html 
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2.2.2. Taeduk Radio Astronomy Observatory Observations 

Observations of CO 1-0 were carried out in 2013 February with the 14 m radio 
telescope of Taeduk Radio Astronomy Observatory (TRAO), equipped with 9 receivers 
arranged in 3 x 3 pattern with 50 " spacing. The central position of target was the same 
as in the KVN observations. The system temperatures during observations ranged from 
950 K to 1100 K depending on the weather conditions and source elevations. The velocity 
resolution is 0.327 km s“^ and the averaged 1 a rms noise level is about 0.37 K. 

The beam FWHM and beam efficiency are 50 " and 0.48 at 115 GHz, respectively. 
We mapped ~275"x 275" region around L1251A, sampling every a half of beam FWHM 
(25"). We also used CLASS reduction package from IRAM for data reduction. A summary 
of observations is also presented in Table 1. 


2.2.3. Data from Other Observations 


The 350 /im continuum emission map fjWu et al.ll2007() of L1251-C was obtained using 
the SHARC-H camera on the Caltech Submillimeter Observatory (CSO). 


The 850 pm continuum emission map was reprocessed by 


Di Francesco et ah 


(1200811 


with the data obtained from the James Clerk Maxwell Telescope (JCMT) Submillimeter 
Common User Bolometer Array (SCUBA). The data used in this work are from the public 
archive maintained by the Canadian Astronomical Data Centrqj (CADC). 


^The Canadian Astronomical Data Centre is operated by the Dominion Astrophysical 


Observatory for the National Research Council of Canada. 











2.3. Radio Interferometer Observations 


2.3.1. Submillimeter Array Observations 


The Submillimeter array (SMA; 


Ho et al. 


2 OO 4 J) observations of L1251A were carried 
out on 2007 October 17 in the compact conhguration with six antennas (6*fwhm ~3 "). 
The quasar J2005+778 and J0102+584 were observed for the gain calibration. Titan was 
observed for the flux calibration, and the quasar 3C 454.3 was observed for the pass-band 
calibration. Natural weighting was used without any tapering to produce the data cubes. 
The 230 GHz receiver was tuned to observe CO 2-1 and ^^CO 2-1 in separate correlator 
windows of 2 GHz bandwidths. Channel widths for CO 2-1 and ^^CO 2-1 are 1.06 km s“^ 
and 0.28 km s“^, respectively. The rms noise levels of CO 2-1 and ^^CO 2-1 are presented 
in Table 2. The remaining chunk of 24 windows were assigned to observe the 1.3 mm 
(225.404 GHz) continuum emission. The continuum image was processed with appropriate 


primary beam correction (the primary beam is ~55 " at this frequency). The 1 a rms 


of 1.3 mm continuum map is 4.5 mJy beam The IDL-b ased MIR so 


adapted for SMA was used for calibration and the MIRIAD fjSanlt et al 


twar e package 


evel 


19951) was used for 


imaging and analysing. A snmmary of observations is shown in Table 2. 


'^See http://cfa-www.harvard.edu/~cqi/mircook.html, 
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3. RESULTS 


3.1. Infrared Images 


3.1.1. Central Triple System: L1251A 


IRAC 3.6, 4.5, 5.8 /xm and MIPS 24 /xm images of L1251A are presented in Fignre 1. 
Three IR sonrces, IRS 1, IRS 2, and IRS 3, were identified with the HDR mode observations 
in the IRAC band images. On the other hand, they are unr esolved in the MIPS band 
image fPignre Id). Mnlti components of L1251A at NIR were detected for the first time by 


Rosvick fc Davidgel (jl995[l , where five NIR sonrces were identified as la 


Deled A-E near IRAS 


22343 +7501. Unfortnnately the coordinates for the five sonrces from 


Rosvick fc Davidge 


(119951) are not snfficiently precise to associate them with the Spitzer sonrces. Comparisons 


of the relative positions show that Spitzer s onrces IRS 1. IRS 2 . and IRS 3 correspond to 


sonrce D, sonrce 4/B/C, and sonrce 1/A of iRosvick fc Davidgel (119951) . respectively. IRS 1 


is the brightest sonrce among three components in the IRAC band images except in the 
IRAC 1 band (3.6 +m), while IRS 2 is the faintest one. 


3.1.2. YSO Candidates in L1251-C 


Fignre 2 shows the IRAC three-co 
For identification of YSO candidates. 


or composit e 


Evans et ah 


imag e of the 3.6, 4.5, and 8.0 pm bands. 


(120071) exclnded possible extra-galactic 


backgronnd sonrces based on the criteria derived with the color-magnitnd e diagrams of 


the r eprocessed Spitzer Wide-area Infrared Extragalactic Survey (SWIRE; 


Lonsdale et ah 


20031) data. According to the c2d YSOc catalog, 16 YSO candidates were identihed in 
this region. Three more sources were identihed by HDR mode observations (IRS 1, IRS 
2, and IRS 3). Includingly added the three sources, 19 YSO candidates were identihed in 
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total, and Figure 2 shows our identification numbers for the sources around L1251-C on 
IRAC 4.5 /rm image. Labels 1,2, and 3 represent IRS 1, IRS 2, and IRS 3, respectively, and 
labels from 4 to 19 are assigned along with Right Ascension. Coordinates and Spitzer flux 
densities of the YSO candidates identihed in the Spitzer images are presented in Table 3. 


3.2. Continuum Images, Molecular Line Spectra, and Maps 


3.2.1. Radio Single-dish Continuum Emission 


Figure 3 shows the continuum images observed with CSO SHARC II at 350 pm and 
JCMT SCUBA at 850 pm. Dust c ontinuum condensations associated with L1251A have 


been detected in both wavelengths (lYoung et ah 


2006 


Wu et ah 


20071 ). With the effective 


beam FWHMs of ~9 " and ~22.9 ", respectively, YSOs in L1251A were not resolved. The 
peak positions of dust continuum emission in both wavelengths are located at the center of 
the three IR sources, suggesting that the dust continuum emission traces the dense core 
covering L1251A with peak intensities of 5.4 and 1.02 Jy beam“^ at 350 and 850 pm, 
respectively. Dust continuum condensations in both wavelengths are seen elongated in the 
east-west (E-W) direction from L1251A. The 850 pm dust continuum emission is likely 
optically thin and traces the dense and cold material within the molecular core. Therefore, 
the total (gas+dust) mass of L1251A, Mtotai, was calculated with the 850 pm integrated 
flux, Fj,, using 


Yftotal — 




and assuming that all dust-continuum emi ssion at 850 nm arises from dust and is optically 


thin. Here, d is the distance (300 ± 50 pc; 


KunetaL 


19931) . and Bi, is the Planck function 


for a blackbody of dust temperature Tjust- To estimate mass, we ass umed Tnnst = 20 K and 


adopted the mass opacity coefficient, kx = 0.1(250pm/A)^ cm^ g ^ flHildebrand 


1983 ) with 
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(3 = 1.8 flHatchell et al 


20131) ■ This mass opacity assumes a gas-to-dust ratio of 100, and 


the estimated mass represents the total (gas+dust) mass of the molecular core. Since 350 
/im is not in the Rayleigh-Jeans limit for this temperature regime, so we derived total mass 
only with the 850 /xm flux. The derived total mass is 4.98 Mq with the integrated flux 
density, 6.18 Jy. The aperture size used for the flux density is 60 ". 


3.2.2. Radio Single-dish Molecular Line Emission 

Figure 4 presents the integrated intensity maps of the HCO’*' 1-0, 1-0, N 2 H”^ 

1-0, and HCN 1-0, which all trace dense gas, made from the KVN single dish observations. 
We integrated over the full widths of HCO’*' 1-0 and 1-0, and all components of 

hyperhne structures were integrated for the integrated intensity maps of N 2 H''' 1-0 and 
HCN 1-0. The molecular line integrated intensity extends in the E-W direction in a manner 
similar to the distribution of dust continuum emission except HCO’*' 1-0, which may trace 
both infall and outflow (Figure 3). The integrated intensity peaks of these molecular lines 
are offset from the continuum peak position of L1251A, and relatively optically thin lines 
such as 1-0 and N 2 H’'' 1-0 show bigger offsets. The integrated intensity peak 

of 1-0 is ~30 " to the south of L1251A (Figure 4b) and is coincident with the 

position of IRS 14. The integrated intensity peak of N 2 H+ 1-0 is ~20 " to the southwest of 
L1251A (Figure 4c). The integrated intensity peaks of relatively optically thick molecular 
emission such as HCO’*' 1-0 and HCN 1-0 slightly shifted to southwest. The disagreements 
of positions between the intensity peaks of single dish dust continuum emission and the 
molecular line integrated intensity peaks might indicate that the distribution of molecular 
gas around L1251A is complicated because of 4 sources existing within about 25 " (i.e., 0.04 
pc). For example, heating by these sources can evaporate CO from the dust grain surfaces, 
which is a mother molecule of HCO"^ but a destroyer of N 2 H+. 
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Figure 5 shows the line prohles of the above molecular lines at the center position (IRS 
1), as well as that of H 2 CO 2i2-lii. The spectrum of HCO’*' 1-0 toward L1251A shows 


a strong self-absorption dip wi th a classic 


believed as an infall indicator (jZhou et ah 


jlue asymmetric line prohle, which has been 


19931) . These dips imply true self-absorption 


when the velocity of the line peak of optically thin tracers such as 1-0, H 2 CO 

2 i 2 - 1 ii, and N 2 H+ 1-0 coincides with the velocity of dip; although the centroid velocities 
of optically thin tracers such as 1-0 and H 2 CO 2 i 2 -lii are slightly blueshifted, 

they peak at the velocity of the self-absorption dip. Gaussian £t FWHMs line widths of 
1-0 and H 2 CO 2i2-lii are 1.3 and 1.9 km s“^, respectively. From simultaneous 
hyperfine structure fitting, the FWHMs of N 2 H+ 1-0 and HCN 1-0 are 0.6 and 1.04 km 
s“^, respectively. Therefore, N 2 H’'' 1-0 seems to trace the most quiescent gas component. 
In addition, the profile of F = 2-1 hyperfine component of HCN 1-0 line, which also has 
been believed as an infall tracer, is similar to that of HCO’*' 1-0. 


The line prohle of HCO’*' 1-0 is superimposed by two (broad and narrow) components 
and the narrow component shows the self-absorption feature. The broad component (ucen ~ 
-4.8 km s“^, Av ~ 2.9 km s“^) is seen in all positions across the observed map while 
the narrow component (ucen ~ -5.3 km s“^, Av ~ 0.8 km s“^) is detected just around 
L1251A. The integrated intensity map of the broad component (not shown) shows a similar 
distribution to the HCO'*' 1-0 outhow (see section 4.4 below), suggesting the origin of the 
broad component may be the outhows from L1251A. The integrated intensity map of the 
HCO'*' 1-0 narrow component and that of the isolated component of N 2 H^ 1-0 show similar 
distributions, although the integrated intensity peaks of the two lines are not coincident 
(Figure 6). The similarity suggests that the narrow component of HCO’*' 1-0 may trace the 
dense and relatively quiescent molecular gas around L1251A like the other molecular lines 
(except for CO 1-0). 
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3.2.3. Radio Single-dish Water Maser Emission 


H 2 O maser emission at 22 GHz was monitored toward L1251A to find variations of 
the systemic velocity and intensity in time. In Figure 7, multi-epoch results are presented. 
Maser components of different velocity or intensity were seen at different epochs. Five 
redshifted maser components can be identified at -2.4, 0, 2, 3.2, and 153 km s“^. The 
-2.4 km s“^ component was presented in 2013 February, slightly weakened in 2013 March 
and become very strong in 2013 May. The 0 km s“^ component was first detected in 2012 
December and became stronger in 2013 January but weakened in 2013 May. The 3.2 km s“^ 
velocity component and the very high velocity component 153 km s“^ were also detected in 
2013 February (Figure 8). The integrated intensity peak of the H 2 O maser emission (not 
shown) located in L1251A but the beam FWHM is not small enough to identify the specific 
location of the water maser source. 


The isotropic luminosity of an H 2 O maser can 


3 e estimated by the integrated intensities 


of observed emission using the following equation flBae et al 


201 ih : 


-^H20 = 


F,.dv. 


( 2 ) 


Here, d is the distance (300 ± 50 pc; 


Kun et al. 


19931 1. and u is the observed frequency. 


When multiple lines are present, the sum of all the individual lines are used to calculate the 
luminosity. The derived H 2 O maser luminosities are presented in Table 4. 


3 . 2 . 4 . Radio Interferometer Continuum Emission 


The 1.3 mm continuum image obtained with the SMA toward L1251A is presented 


in Figure 9. Three individual condensation s were detected. For 


sources, we will follow the nomenclature of 


Reipurth et al. 


he 1.3 mm continuum 


f 2004j l. Although the peak 


positions of submillimeter sources are not exactly coincide with the positions of IR sources. 
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VLA 6 and VLA 7 are the same sources as IRS 1 and IRS 3, respectively. However, VLA 
10 is apart from IRS 2 by about 2 ". The positions of IRS 2 and the 2MASS source are 
consistent. Therefore, VLA 10 and IRS 2 are tracing two different s ources. The du st 
continuum emission related to VLA 6 and VLA 7 was first detected by 
their VLA A and VL A B) at 3.6 cm and next by 


Meehan et ah 


Beltran et af 


ill998 


Reipurth et ah 


fl200ll ) at 3.6 cm and 6 cm. 


(12004[) . however, first detected the continuum condensation associated with 


VLA 10 at 3.6 cm, and they mentioned that VLA 10 is a radio jet source at a position angle 
of about 53 °, indicating that VLA 10 is anoth er protostellar membe r in L1251A although 


no corresponding IR source has been detected. 


Reipurth et ah 


(120041) also mentioned that 


VLA 7 is a radio jet source at a position angle of about 126 °. 

Although VLA 6 and VLA 10 are connected at low intensity levels, three sources are 


resolved in our 


.3 mm continuum map. To estimate their masses, we used T^ust = 30 K 


(Joreensen et ah 

2007 

) and ac 

optec 

. the dust mass opacity coefficient, kx = 0.02(1 mm/X)^ 

), which is normalized for flat sources (B = 1). The 

cm^ g ^ ( 

Beckwith & Sarerent 

1991 

masses of VLA 6, VLA 7, and VLA 10 are presented in Table 5. The total mass derived 


from the SCUBA 850 /im continuum map is roughly 55 times larger than the sum of masses 
derived from the SMA 1.3 mm continuum sources because the single dish flux includes 
extended emission around all four sources. 


3.2.5. Radio Interferometer Moleeular Line Emission 

Line profiles of CO 2-1 and ^^CO 2-1, observed with the SMA at the position of each 
submillimeter source are presented in Figure 10. Channel widths for CO 2-1 and ^^CO 2-1 
are 1.09 km s“^ and 0.23 km s“^, respectively. Among the three submillimeter sources, the 
strongest peak brightness temperature of CO 2-1 (24.75 K) is seen at VLA 6, while the 
peak brightness temperatures of ^^CO 2-1 (~7 K) are similar toward all three sources. The 
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centroid velocity of 2-1 is shifted from red to blue in the order of VLA 6, VLA 10, 
and VLA 7. There are channels without emission in 2-1 at uisr = -4.94 and -3.84 km 
s“^. The two dips of 2-1 cannot be resolved with the ~1 km s“^ velocity resolution of 
CO 2-1 data. Moreover, CO 2-1 emission is missing in the velocity range of wisr = -6 to -4 
km s“^. At these velocities, most emission from the larger-scale structure in that velocity 
range has been likely resolved out. 

In Figure 11, the velocity channel map of the CO 2-1 line is presented. The CO 2-1 
emission is in the velocity range of -14.9 to 5.2 km s“^. In the channel map, the CO 2-1 
emission is missing at the channels of -5.4 and -4.3 km s“^. Moreover, the strong negative 
contours are shown near bright emission, caused by the missing flux effect at -3.3 and -2.2 
km s“^. A blueshifted component related to VLA 6 or VLA 10 is seen extending to the 
southwest of L1251A. A strong redshifted component is located at the position of VLA 6, 
extending in the northeast direction. Though it is hard to distinguish apparent bipolar 
outflow structures, we infer that the blue- and redshifted components are from outflowing 
gas, because the velocity offsets from the systemic velocity are > 2.5 km s“^. We will 
discuss more about outflows in section 4.4 below. Figure 12 shows the CO 1-0 integrated 
intensity across L1251A, overlaid onto the 1.3 mm continuum emission map for reference. 
The integrated intensity map of CO 1-0 shows E-W elongation, and the intensity peak is 
located at the position of VLA 6. 

Figure 13 shows the velocity channel map of the ^^CO 2-1 line. Every four channels 
are integrated for ^^CO 2-1 to compare with the CO 2-1 channel map. The redshifted 
emission at the position of VLA 6 shown at -3.1 km s“^, while the blueshifted emission at 
the position of VLA 10 and VLA 7 is shown at -5.4 and -6.5 km s“^, respectively. The 
blueshifted emission at -5.4 km s“^ is slightly extended to the southwest and the blueshifted 
emission at -6.5 km s“^ is extended to the northeast. At -5.4 km s“^, there is another blob 
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to the east of VLA 10, which can be discerned in the integrated intensity map of 2-1. 
In Fignre 14, the integrated intensity map of 2-1 across L1251A is presented and also 
overlaid onto 1.3 mm continnum emission map for reference. Each continunm sonrce seems 
to be associated with a 2-1 integrated intensity peak althongh those peaks are not 

completely distinct. The whole emission structnre extended along the NE-SW direction. 

Most emission of blue- and redshifted components of 2-1 is seen in the velocity 
range where the CO 2-1 emission appears to have been resolved out. In addition, the 
blue- and redshifted components of 2-1 have low velocities (within uisr ± 3 km s“^) 
compared to that of CO 2-1 (within uisr ± 9 km s“^). Therefore, the features of ^^CO 2-1 
seem to trace the innermost part of envelopes around each submillimeter source rather than 
outflows. 


4. ANALYSIS 


4.1. YSO Classification 


Except IRS 1, IRS 2, and IRS 3 identihed by the HDR observations, the IR spectral 
index a of 16 YSO candidates are listed in the c2d YSOc catalog. The a for IRS 1, IRS 2, 
and IRS 3 are calculated in this work. All photometries b etween 2 nm and 24 pm in the 


c2d catalog are used to determine the IR spectral index a flEvans et ah 

d log(A^(A)) 


200911 : 


a = 


Here, A is the waveleng 


criteria of 


Evans et ah 


d log(A) 

h ari d S'(A) is the flux density at that wavelength. Following the 


( 3 ) 


f 20091 1 (i.e.. Class I for 0.3 ^ a. Flat -0.3 ^ a < 0.3, Class II for 


-1.6 ^ a < -0.3, and Class III for a < -1.6), IRS 1, 11, and 16 are Class I candidates, IRS 
3, 7, and 17 is Flat class candidates, and all others are Class II candidates. Table 6 lists 
the IR spectral index determined for each source. We calculated the IR spectral indices 
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of the newly added YSOs, also listed in Table 6. We also applied two other criteria of 
classihcation (the color-color diagram and Tboi) to all 19 YSOs. Photometric data of all 
YSOs are presented in Table 3 from the Spitzer observations. Figure 15 shows a [3.6] - 
[4.5] versus [5.6] - [8.0] color-color diagram used to classify the evolutionary stages of YSO 
ca ndidates, includi ng all identihed sources in Figure 2. We ado pted the same crit eria used 


m 


Lee et ah 


(120061 ). which is basically based on the models of 


Allen et ah 


(120041) . Two 


sources, IRS 1 and 14, are identihed as Class O/I candidates in color-color diagram. IRS 11 
and 17 are located at the color-color boundaries between Class O/I and Class Ih The IR 
spectral indices of these sources, however, indicate that IRS 11 is in Class I and IRS 17 is 
in Flat class. IRS 4 is located close to the boundary of Class III criteria, but its IR spectral 
index of IRS 4 is in the range for Class Ih 

and bolometric temperatures (Tboi) 


In Table 6, derived bolomet ric luminosities (Lhoi 

' 3 


are listed. If the classihcation by 


Mvers fc LaddI (119931) (i.e.. Class 0 for Tboi < 70 K, Class I 


for 70 K ^ Tboi ^ 650 K, and Class II for 650 K < Tboi ^ 2800 K) is used, IRS 1, 11, and 16 
are identihed as Class I, and the other 16 YSOs are identihed as Class Ih The evolutionary 
stages of most sources identihed by Tboi agrees with those by the IRAC color-color diagram. 
Note, however, that IRS 14 is identihed as an embedded source by the color-color diagram 
but as Class II by Tboi and a. It is more likely in Class II, but since it is spatially located 
at the H^^CO’*' 1-0 emission peak it may be Class O/I instead. Furthermore, IRS 11, which 
is located at the boundaries of Class O/I and Class II in the IRAC color-color diagram, 
is identihed as Class I by Tboi (IRS 11 has a very low luminosity). In addition, IRS 4 is 
identihed as Class II source by Tboi. Finally, IRS 16 is identihed as a Class I source with 
the lowest Tboi among 19 YSO candidates. 


IRS 11 and IRS 16 are identihed as Class I YSOs with luminosities smaller than 0.1 
Lq, suggesting these sources are candidates of Very Low Luminosity Objects (VeLLOs, e.q.. 













18 


Youne; et al. 

to 

o 

o 

_ 

Bourke et al. 

2005 


Dunhametal 


20061) . Ho wever, it cannot be ruled 


out that these two sources may be background galactic sources. 


Evans et ah 


( 20071) dehned 


the probability of being background galactic source, Log(Prob_Galc), and galaxy candidates 
have Log(Prob_Galc) ^ -1.47. For IRS 11 and IRS 16, the probability is -2.17 and -1.58, 
respectively. 


Spectral energy distributions (SEDs) for the YSOs in L1251-G are presented in Figure 
16. The SEDs include 2MASS, IRAG, and MIPS bands data. IRS 1 is the brightest source 
at near- and mid-infrared wavelengths except in IRAG 1 band (3.6 pm). 


4.2. L1251A and the Common Envelope 


Including the radio continuum source VLA 10, the number of members in L1251A 
is four. The mean projected separation between each companion is about 1220 AU. 


Looney et ah 


fl2000f) classihed multiple systems by their morphologies into three groups: 


separate envelope system, common envelope system, and common disk system. The typical 
separation of common envelope systems is from 150 AU to 3000 AU. The projected distances 
among members of this system suggest that L1251A is a common envelope system, having 
one primary core of gravitational concentration. 


4.3. Molecular Outflows 


In L125 1-G, an exten d ed bip olar GO 1-0 outflow driven by IR AS 22343-1-7501 was 


reported by 


Sato fc Fukuil (119891) . An optical jet was also observed (iBalazs et al 


19921) 


with the same axis ( NE-SW) as the extended GO 1-0 outflow. The outflow map of 


Schwartz et al. 


(119881) . however, shows a bipolar structure with an unclear outflow axis in a 


small region around L1251 A. 
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For this paper, the velocity ranges for blue- and redshifted outflow wings were dehned 
as follows. The outer velocities of the ranges are determined for each position where the 
emission is greater than the respective la rms noise level. Spectra at outflow-free positions, 
chosen by eye, are used to dehne the inner velocities of the ranges, as being equal to the 
outermost velocities. 


The averaged outflow-free spectra of CO 1-0 for the blue and red component are 
presented in Figure 17a. The spectra at the peak positions of the blue- and red-outflow 
components are also shown. The outflow map of CO 1-0 made by integrating intensity 
over the determined outflow wings is presented in Figure 17b. The outflow wing emission 
is elongated mainly in the east-west (E-W) direction, and it also shows asymmetry in 
intensity and morphology. For example, the red component of the CO 1-0 outflow is more 
broadly extended than the blue component. Also, the blue peak is located at the east of 
L1251A, and the red peak is located at the north of that region. 


In comparison, a HCO’*' 1-0 outflow map shows blue and red components with 
signihcant asymmetry in intensity and size. The velocity ranges for blue- and redshifted 
outflow wings of HCO^ 1-0 are -9.7 to -5.9 km s“^ and -4.0 to 0.2 km s“^, respectively. 
Figure 18 shows the blue component is dominant. A poorly collimated outflow of HCO^ 
1-0 is extended along the NW-SE direction, i.e., inconsistent with the E-W direction of 
the CO 1-0 outflow. The position angle of the HCO+ 1-0 outflo w axis, however, doe s 


correspond to that of the HCO’*' 1-0 outflow-like structure seen by 


Nikolic et ah 


fl2003h . 


Figure 19 shows the outflow of CO 2-1 detected with the SMA. The velocity ranges 
for blue- and redshifted outflow wings of CO 2-1 are -20 to -7 km s“^ and from -2 to 
10 km s“^, respectively. Then the axis of this small scale outflow is not consistent with 
either the CO 1-0 outflow or the HCO^ 1-0 outflow. The blue- and redshifted emission 
extends along the NE-SW direction and is very compact. IRS 1 located at the position of 
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the red peak of the CO 2-1 outflow and IRS 1 is the closest to the center of the blue and 
red components. In addition, IRS 1 is in the earliest evolutionary stage among three YSOs 
of L1251A detected in IRAC observations (see 5.1 below), suggesting IRS 1 may be the 
driving source of the compact CO 2-1 outflow. However, in spite of no corresponding IR 
source, a radio jet has been detected in VLA 10 and its propagation direction is consistent 
with the compact outflow direction. Therefore, higher-resolution observations are needed to 
determine the engine of this compact outflow. 


Different molecular outflows on different scales show dissimilarities in axis, morphology. 


and strength (Figure 17, 18, and 19), sugges 
L1251A. The extended CO 1-0 outflow from 


ing that ther e are m ultiple outflows from 


Sato fc Fukuil fjl989l. see there Figure 3) has 


a direction completely opposite to that of the compact CO 2-1 outflow presented here. 
The HCO’*' 1-0 and CO 1-0 outflows presented here may show a mixture of extended and 
compact outflows with poor collimation. 


5. Discussion 

5.1. Star Formation in Small Groups 

In L1251-C, 19 YSOs are detected, including three that are identihed in the HDR mode 
observations. The submillimeter source, VLA 10, is excluded for this analysis, because VLA 
10 was not detected by Spitzer observations. The fraction of identihed sources as Class O/I 
YSOs is 16 % and 84 % of sources are identihed as Class II, based on Tboi classihcation 
(see Table 6). In L1251A, however, one is a Class O/I source and two others are Class II 
sources. The whole held size of the MIPS observations of L1251-C is ~1.44 pc^, therefore 
the surface density of YSOs in L1251-C is ~13.2 pc“^. The surface density of YSOs in 
L1251A is larger than that of the whole L1251-C region by two orders of magnitude. The 
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area used for calculating the YSO surface density of L1251A is ~0.08 pc^, which covers its 
emission at submillimeter wavelengths (see Figure 3). 


We comp are L1251-C wit h the other core in L1251 (jSato et ah 


with Spitzer (ILee et ah 


20061). 


Lee et ah 


199411 and also observed 


fj2006l) shows that the surface density of YSOs 


in L1251-E is close to those in active star-forming regions. L1251-E (~2.5 pc^) has a 
larger area than L1251-C (~1.44 pc^). In L1251-E, 20 YSO candidates (4 Class O/I, 14 
Class II, and 3 unclassihed YSOs) were detected in the MIPS 24 /im band. As a result, 
in L1251-E, the surface density of YSOs is ~10 pc“^, which is slightly lower than that of 
L1251-C by a factor of 1.3. L1251B, the densest part in L1251-E, however, shows a larger 
YSO surface density than that of L1251A by a factor of 2. The composition of YSOs in 
various evolutionary stages in L1251-E and L1251-C shows similar fractions. For example, 
the fractions of Class O/I candidates in L1251-C and L1251-E are about 16 % and 22 %, 
respectively. The fraction of Class O/I sources in L1251B (50 %), is slightly larger than 
that in L1251A (~34 %). Overall, L1251-C has a similar star formation environment to 
L1251-E. 


One Class O/I YSO and two Class II YSOs are located in L1251A in the IRAC HDR 
images (see, e.g.. Fig. 1). The average projected distance among the 19 YSOs in L1251-C 
is about 136 ", while the average projected distance of three YSOs in L1251A is about 
5.6 ", suggesting clustering. The intensity peak of single-dish submillimeter continuum 
emission, which is coincident with all three YSOs of L1251A, is located between all three 
YSOs, suggesting that three YSOs at different evolutionary stages all reside in a singe 
dense core. Interferometer submillimeter molecular line emission from L1251A observed 
with SMA, which shows more detailed view of the dense gas around the YSO group and 
resolves out the large r-scale structure, has their intensity peaks located at the position 
of IRS 1, showing IRS 1 is more deeply embedded than the other two YSOs. Single-dish 
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submillimeter molecular line emission from L1251A, however, has their integrated intensity 
peaks shifted to the south of L1251A, probably because of CO evaporation. 

The interferometric continuum data at 1.3 mm are likely tracing disks around 
protostars, but the contamination from the inner envelope cannot be ruled out because 
the deconvolved sizes of submillimeter sources are similar or slightly larger than the beam 
FWHM (2.7 "x 2.5 ”). The masses derived from the 1.3 mm continuum decrease in order 
of VLA 10, VLA 7, and VLA 6 (0.030, 0.023, and 0.020 M©, respectively). VLA 10 does 
not have a corresponding IR source, indicative of a deeply embedded source. Therefore, the 
1.3 mm continuum emission of VLA 10 might trace its very dense inner envelope. 


5.2. Outflows in L1251A 


At least two outflows with d ifferent properties exist in L1251-C. One is an extended CO 


1-0 outflow f Sato fc Fukui 


The extended outflow from 


989h and the o ther is a compact CO 2-1 outflow (this paper). 


Sato fc Fukuil (119891. see their Figure 3) shows a completely 


opposite outflow direction to that of the compact outflow from this work. The dynamical 
time of the extended outflow (1.7x10^ yr) is greater than that of the compact outflow 
(4.7x10^ yr) by two orders of magnitude. Submillimeter continuum associated with radio 
jets (VLA 7 and VLA 10) is detected in L1251A. VLA 7 has the NW-SE propagation 
direction and VLA 10 has the NE-SW propagation direction. Poorly collimated distribution 
of blue- and redshifted components of KVN HCO’*' 1-0 outflow and TRAO CO 1-0 outflow 
may be also affected by these radio jets of VLA 7 and VLA 10. 


The driving sources of each outflow, however, are hard to distinguish. A strong 
candidate of the driving source of the compact outflow is IRS 1 (VLA 6) because it is an 
embedded Class I source, and the blue- and redshifted lobes are closely associated with the 
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position of IRS 1 (see Figure 19). However, VLA 10 cannot be ruled out because its dust 
continuum emission seems to trace a very dense inner envelope (indicative of an extremely 
embedded source), and it is associated with a radio jet, whose p ropagation direction is 


consistent with that of the compact outflow flReipurth et ah 


20041 . 


6. SUMMARY 

L1251-C was observed as part of the Spitzer legacy project, c2d. Also, various 
molecular emission lines associated with L1251A have been observed with the KVN 
21 m radio telescope and TRAO 14 m radio telescope to understand its environment. 
Single-dish continuum data from other observations are also presented. In addition, SMA 
interferometric continuum and line data of L1251A were obtained at high angular resolution. 
With these data, we hud: 

1. 19 YSO candidates are associated with L1251-C, including three additionally 
identified in the HDR mode Spitzer IRAC bands: IRS 1, IRS 2, and IRS 3. 

2. Single-dish continuum emission at 350 and 850 /xm is extended to the E-W and may 

trace an envelope containing all three sources. Individual condensations associated with 
each IR source are not resolved in the dust continuum map at 350 and 850 /xm. L1251A 
was mapped at KVN in HCO'*' 1-0, 1-0, N 2 H''' 1-0, and HCN 1-0. The integrated 

intensity peaks of all molecular lines observed at KVN are offset from the peak positions 
of the continuum maps at 350 and 850 /xm. Notably, the integrated intensity maxima of 

1-0 and N 2 H''' 1-0, which are optically thin, are offset by tens of arcseconds from 
L1251A to the south or southwest. The line profile of a well known outflow and infall tracer 
HCO'*' 1-0 line are composed of two velocity (broad and narrow) components and show a 
self-absorption feature. Though this feature may trace infall, the broad component is a 
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sign of outflow and the narrow component may trace dense envelope material (like other 
molecular lines described here). Several velocity components of the H 2 O 22 GHz maser 
have been detected during multi epoch observations. 

3. Three dust continuum condensations (VLA 6, VLA 7, and VLA 10) are detected 
by the SMA observations at 1.3 mm. Two of three (VLA 6 and VLA 7) are detected at 
the positions of IR sources (IRS 1 and IRS3, respectively) in L1251A. The positions of IRS 
2 and VLA 10 are apart by ~2 ", suggesting they are different sources. Finally, the total 
number of protostars in L1251A is four by adding the most embedded source of VLA 10, 
which does not have a corresponding IR source. The SMA integrated intensity maps of CO 
2-1 and 2-1 show complicated distributions although a compact outflow is clearly 

detected in CO 2-1 in the NE-SW direction. 


4. In L1251-C, IRS 1, IRS 11, and IRS 16 are identihed as Class I, and the other 16 
YSOs are identihed as Class 11. The surface density of L1251-C (~1.44 pc^) is similar to 
that of L1251-E (~2.5 pc^). The fraction of Class O/I sources (16 %) in L1251-C is also 
similar to that seen in L1251-E (22 %). 


5. Molecular outhows on different scales have been mapped in HCO’*' 1-0, CO 
1-0, and CO 2-1. These outhows lack coincidence in axis, morphology, and strength, 
sugg esting there ar e two outhows from two sources. One is a extended CO 1-0 outhow 


from 


Sato fc Fukuil fjlOSOl) and the other is a compact SMA CO 2-1 outhow described 


here. KVN HCO"*" 1-0 and TRAO CO 1-0 observations may show a mixture of extended 
and compact outhows. IRS 1 (VLA 6) is likely the driving source of the compact outhow 
detected in higher-resolution SMA CO 2-1 observations, but VLA 10 cannot be ruled out. 
Therefore, interferometric observations with still higher spatial resolution are necessary to 
identify better which source is driving the compact outhow. 
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Table 1. Summary of Radio Single-Dish Observations 


Molecular line 

V 

[GHz] 

A V 

[km s“^] 

Beam FWHM 

[arcsec] 

'f]mb 

1 a rms 

[K] 

Telescope 

H 2 O 016-523 

22.235080 

0.21 

120 

0.38 

0.08 

KVN (YS) 





0.42 

0.09 

KVN (US) 

H13CO+ 1-0 

86.754288 

0.22 

31.5 

0.42 

0.08 

KVN 

HCN 1-0 

88.631847 

0.22 

31.5 

0.41 

0.09 

KVN 

HCO+ 1-0 

89.188526 

0.22 

31.5 

0.40 

0.09 

KVN 

N 2 H+ 1-0 

93.176265 

0.22 

31.5 

0.37 

0.09 

KVN 

H 2 CO 2i2—111 

140.839520 

0.27 

23.4 

0.30 

0.11 

KVN 

CO 1-0 

115.271202 

0.21 

50 

0.48 

0.37 

TRAO 
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Table 2. Summary of SMA Observations 


Molecular line 

V 

[GHz] 

Side Band 

Number of 

Channels 

A V 

[km s“^] 

Beam FWHM 

[arcsec x arcsec] 

P.A. 

[deg] 

1 cr rms 

[Jy beam“^] 

CO 2-1 

230.5379700 

USB 

128 

1.06 

3.6 X 3.1 

-1.3 

0.20 

13CO 2-1 

220.3986765 

LSB 

512 

0.28 

3.7 X 2.3 

-5.8 

0.37 






Table 3. Infrared Fluxes of Young Stellar Object Candidates from the Spitzer Data 


Source 

RA 

jh m sj 

Dec 

[O , 

IRAC 3.6 //m 

[mjy] 

IRAC 4.5 fxm 

[mJy] 

IRAC 5.8 fira 

[mJy] 

IRAC 8.0 /J.m 

[mJy] 

MIPS 24 fim 

[mJy] 

MIPS 70 fim 

[mJy] 

1 

22:35:23.28 

+75:17:07.73 

293.00 (13.00) 

748.00 (13.00) 

1,550.00 (32.00) 

2,510.00 (39.900) 



2 

22:35:24.57 

+75:17:03.61 

183.00 (3.20) 

114.00 (3.60) 

111.00 (3.30) 

128.00 (1.30) 



3 

22:35:24.96 

+75:17:11.33 

532.00 (22.00) 

610.00 (11.00) 

707.00 (14.00) 

890.00 (19.00) 



4 

22:34:12.00 

+75:18:09.36 

183.00 (58.90) 

114.00 (14.60) 

111.00 (34.90) 

128.00 (31.40) 

1,440.00 (135.00) 

1,740.00 (162.00) 

5 

22:34:40.56 

+ 75:17:44.16 

30.20 (1.51) 

28.40 (1.38) 

27.70 (1.36) 

35.60 (1.40) 

139.00 (12.90) 

257.00 (24.50) 

6 

22:34:48.24 

+ 75:13:35.04 

8.97 (0.44) 

7.30 (0.35) 

6.68 (0.32) 

9.04 (0.43) 

17.80 (1.66) 

50.90 (5.83) 

7 

22:35:00.72 

+75:15:36.36 

22.00 (1.13) 

23.40 (1.13) 

25.70 (1.25) 

32.30 (1.52) 

102.00 (9.44) 

146.00 (14.50) 

8 

22:35:02.40 

+75:17:58.20 

35.60 (1.83) 

43.10 (2.13) 

45.60 (2.23) 

47.60 (2.30) 

65.30 (6.06) 


9 

22:35:04.08 

+75:18:20.52 

32.40 (1.64) 

32.00 (1.54) 

36.00 (1.74) 

40.80 (1.95) 

80.30 (7.44) 

106.00 (10.70) 

10 

22:35:09.60 

+75:16:09.48 

7.96 (0.40) 

8.40 (0.41) 

6.19 (0.32) 

5.13 (0.24) 

12.70 (1.19) 


11 

22:35:14.06 

+75:15:02.39 

1.50 (0.08) 

1.69 (0.11) 

2.56 (0.12) 

20.00 (1.86) 

146.00 (14.30) 


12 

22:35:16.56 

+75:18:46.80 

75.90 (4.02) 

74.20 (3.81) 

73.00 (3.64) 

85.90 (4.15) 

107.00 (9.96) 

58.50 (6.78) 

13 

22:35:25.44 

+75:17:56.04 

102.00 (5.39) 

108.00 (5.38) 

121.00 (6.05) 

123.00 (5.97) 

75.50 (7.45) 


14 

22:35:26.62 

+75:16:36.86 

34.50 (1.78) 

53.20 (2.64) 

58.00 (2.81) 

54.10 (2.71) 

67.90 (8.19) 


15 

22:35:27.12 

+75:18:01.80 

37.20 (1.85) 

34.20 (1.66) 

33.00 (1.64) 

3.55 (1.70) 

52.90 (5.25) 


16 

22:35:59.21 

+ 75:17:50.40 

0.071 (0.007) 

0.078 (0.007) 

0.084 (0.024) 

0.118 (0.020) 

7.820 (1.760) 


17 

22:36:05.52 

+ 75:18:32.04 

78.20 (4.11) 

102.00 (5.30) 

145.00 (7.29) 

246.00 (12.20) 

534.00 (45.60) 

377. (36.60) 

18 

22:36:35.52 

+ 75:21:34.92 

54.80 (2.72) 

54.90 (2.77) 

58.10 (2.81) 

72.60 (3.65) 

206.00 (19.10) 


19 

22:36:59.04 

+75:21:20.52 

2.87 (0.19) 

2.38 (0.14) 

2.11 (0.13) 

2.06 (0.11) 

2.56 (0.30) 
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Table 4. The isotropic luminosity of H 2 O maser 


Observing Date 

-bUaO 

[Le] 

2012 Dec 

9.46E-09 

2013 Jan 

5.55E-09 

2013 Feb 

7.41E-08 

2013 Mar 

9.59E-08 

2013 May 

2.01E-07 






Table 5. Parameters of the SMA 1.3 mm Dust Continuum Condensations Associated with L1251A 


ID 

RA 

m sj 

Dec 

[O , 

Peak Flux 

[mjy beam“^] 

Integrated Flux Density 

[mJy] 

Source Size 

[" X "] (P.A. ° ) 

Mass 

[Mq ] 

VLA 6 (IRS I) 

22:35:23.45 

+75:17:07.73 

24.20 (0.001) 

35.23 

3.9 X 2.5 (-0.9) 

0.024 

VLA 7 (IRS 3) 

22:35:24.98 

+75:17:12.17 

31.83 (0.003) 

40.53 

3.3 X 2.7 (55.5) 

0.028 

VLA 10 

22:35:24.11 

+75:17:04.81 

38.36 (0.004) 

53.12 

3.3 X 2.8 (-30.6) 

0.037 


oo 

oo 
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Table 6. Bolometric Luminosity, Bolometric Temperature, and IR Spectral Indices of 

YSOs in L1251-C 


Sources 

-^bol 

[^ 0 ] 

Tbol (Class) Spectral index a (Class) 

[K] 

1 

3.2 

622 (I) 

2.94 (I) 

2 

0.54 

1180 (II) 

-0.72 (II) 

3 

2.4 

1041 (II) 

-0.12 (Flat) 

4 

3.1 

1074 (II) 

-0.51 (II) 

5 

0.25 

813 (II) 

-0.38 (II) 

6 

0.085 

1440 (II) 

-0.84 (II) 

7 

0.17 

747 (II) 

-0.29 (Flat) 

8 

0.18 

910 (II) 

-0.59 (II) 

9 

0.12 

960 (II) 

-0.54 (II) 

10 

0.042 

1209 (II) 

-1.02 (II) 

11 

0.015 

367 (I) 

0.35 (I) 

12 

0.47 

1321 (II) 

-0.84 (II) 

13 

0.63 

1422 (II) 

-1.05 (II) 

14 

0.16 

713 (II) 

-0.34 (II) 

15 

0.24 

1441 (II) 

-0.91 (II) 

16 

0.004 

174 (I) 

1.20 (I) 

17 

0.80 

763 (II) 

0.10 (Flat) 

18 

0.35 

1505 (II) 

-0.49 (II) 

19 

0.014 

1537 (II) 

-1.13 (II) 
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Fig. 1.— The Spitzer images of L1251A in log scale at (a) IRAC 3.6 fim, (b) IRAC 4.5 fim, 
(c) IRAC 5.8 yum, and (d) MIPS 24 pm. In Figure Id, Xs indicate the positions of the IR 
sources in L1251A. 
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Fig. 2.— The Spitzer three-color composite image of L1251-C. The IRAC 3.6, 4.5, and 8.0 
/im are presented in bine, green, and red, respectively. Nnmbers (4 to 19) are assigned along 
with Right Ascension after the nnmbers (1 to 3) of three sonrces in L1251A, which are also 
assigned along with Right Ascension. 
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Fig. 3.— Dust continuum emission map t oward L1251A. (a) Continuum emission at 350 


pm obtained with SHARC-II at the CSO f Wu et ah 


20071) . Contour levels are 20, 40, 60, 


80, and 100 % of the peak intensity 5.40 Jv beam fb 

Continuum emission at 850 um 

obtained with SCUBA at the JCMT ( 

Di Francesco et ah 

2008 

). Contour levels are 20, 40, 


60, 80, and 100 % of the peak intensity, 1.20 Jy beam The respective beam FWHM is 
presented in each panner with hlled circle. YSO candidates are denoted by hlled diamonds. 
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RA Offset (arcsec) 


Fig. 4.— Integrated intensity maps of (a) HCO+ 1-0, (b) 1-0, (c) N 2 H+ 1-0, and 

(d) HCN 1-0 on top of the 850 pm map (gray scale). Contour intervals are 10 % of the each 
peak intensity and range from 35 % to 95 %. The gray scale levels are 20, 40, 60, 80, and 100 
% of the peak intensity, 1.20 Jy beam“^. YSO candidates are denoted by hlled diamonds. 
The hlled circle at the bottom-right of each panel denotes the respective beam FWHM. 
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Fig. 5.— Spectra of the all molecular transitions observed with KVN at the center position 
of single-dish molecular line map. The horizontal dashed lines and the vertical dashed lines 
represent the zero baseline and the systemic velocity (-5 km s“^) of L1251-C, respectively. 
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(a) HCC^ 1 -0 narrow component (b) N 2 H'^ 1 -0 isolated component 



RA Offset (orcsec) 

Fig. 6.— Integrated intensity maps of (a) -5.3 km s“^ (narrow) component of HCO’*' 1-0 
and (b) -5 km s“^ (isolated) component of N 2 H+ 1-0. YSO candidates are denoted by filled 
diamonds. Contour intervals are 10 % of the each peak intensity and range from 35 % to 95 
% . The hlled circle at the bottom-right of each panel denotes the respective beam FWHM. 
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-15 -10 -5 0 5 10 15 

V|sr (km s~') 


Fig. 7.— Spectra of the 22 GHz H 2 O maser line toward L1251A. The observation dates are 
labeled. The horizontal dotted line and the vertical dotted line represent the zero baseline 
and the systemic velocity (-5 km s“^) of L1251-C, respectively. 
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Fig. 8.— The highest velocity component of 22 GHz H 2 O maser line toward L1251A. This 
high velocity component was detected only on 2013 Feb. The horizontal dotted line and the 
vertical dotted line represent the zero baseline and the central velocity of the highest velocity 
component, respectively. 
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Fig. 9.— 1.3 mm Dust continuum emission map toward L1251A obtained with SMA. Con¬ 
tour levels are 40, 50, 60, 70, 80, and 90 % of the peak intensity, 0.038 Jy beam“^ and 
the wedge indicates the continuum intensity scale in units of Jy beam“^. The synthesized 
beam FWHM is presented at the bottom-right corner. The rms noise level 1 a is 4.7 mJy 
beam~^. The positions of IR sources in L1251A and the position of VLA 3.6 cm continuum 


sources (IReipurth et ah 


20041 ) are denoted by Xs and triangles, respectively. The size of this 


presented map is the same as the primary beam size at this frequency (~55 
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Fig. 10.— Spectra of CO 2-1 (black) and 2-1 (red) observed with SMA toward 

the three SMA 1.3 mm continuum sources, VLA 6, VLA 10 and VLA 7, from the top to 
bottom. The black dotted vertical lines indicate the systemic velocity of L1251-C, -5.0 km 
s“^, obtained from single dish observations. 
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Fig. 11.— Channel map of CO 2-1 observed with the SMA. Respective velocities in km s“^ 
are indicated in the upper left of each panel. The contours start from 5 a and increase in 
steps of 5 a, where 1 a is 0.268 Jy beam“^. Dashed contours denote negative intensities 
with same steps of positive contours. Crosses denote the SMA 1.3 mm continuum source 
positions (VLA 6, VLA 7, and VLA 10). The synthesized beam FWHM is shown with the 
black ellipse. 




RA offset (arcsec) 

Fig. 12.— Integrated intensity map of CO 2-1 overlaid on the 1.3 mm continuum emission 
(gray scale). The contours represent 30, 45, 60, 75, 80, and 95 % of the peak intensity of 
57.3 Jy beam“^ km s“^ and the wedge indicates the continuum emission scale from 0.015 to 
0.040 Jy beam“^ . The synthesized beam FWHM is shown with the black ellipse. The map 
covers the primary beam size at this frequency (~55 "). 
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Fig. 13.— Channel map of 2-1 observed with the SMA. Respective velocities in km 
s“^ are indicated in the upper left of each panel. The contours start from 2 a and increase in 
steps of 1 a, where 1 a = 0.404 Jy beam“^. Dashed contours show the negative intensities 
with same steps of positive contours. The synthesized beam FWHM is shown with the black 
ellipse. 
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Fig. 14.— Integrated intensity map of 2-1 overlaid onto the 1.3 mm continuum (gray 
scale). The contours represent 25, 40, 55, 70, 85, and 100 % of the peak intensity of 7.75 Jy 
beam“^ km s“^ and the wedge indicates the continuum emission scale from 0.015 to 0.040 
Jy beam“^. The synthesized beam FWHM is shown with the black ellipse. 
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Fig. 15.— [3.6] - [4.5] versus [5.6] - [8.0] color-color diagram for L1251-C YSO candidates in 
the IRAC bands. Source numbers are labeled. The dashed lin e s des cribe the approximate 


criteri a for Class I, Class II, and Class III sources by 


Lee et ah 


(1200411 . 


f 200611 based on 


Allen et ah 
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Fig. 16.— SED of the YSOs in L1251-C. The 
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Fig. 17.— (a) The blue averaged spectrum of the blue wing free positions overlaid on the 
spectrum at the blueshifted emission peak position {top) and the red averaged spectrum 
of the red wing free positions overlaid on the spectrum at the redshifted emission peak 
position (bottom), (b) Integrated intensities of blue- and redshifted emission of CO 1-0. Blue 
contours represent intensities integrated over the blue-shifted range from systemic velocity. 
Red contours show intensities integrated over the red-shifted range. The velocity range used 
for the blue and red components are from -10.4 to -7.1 km s“^ and from -2.2 to -0.2 km 
s“^, respectively. The contours start from 40 % of the peak intensity to 90 % and increase in 
steps of 10 %. IRS 1 and the other YSOs are denoted by a cross and diamonds, respectively. 
The beam FWHM is presented with gray filled circle. 
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Fig. 18.— Integrated intensities of blue- and redshifted emission of HCO’*' 1-0. Blue con¬ 
tours represent intensities integrated over the blueshifted range from systemic velocity. Red 
contours show intensities integrated over the redshifted range. The velocity range used for 
the blue and red components are from -9.7 to -5.9 km s“^ and from -4.0 to -0.2 km s“^, 
respectively. The contours start from 40 % of the peak intensity to 90 % and increase in 
steps of 10 %. IRS 1 and the other YSOs are denoted by a cross and diamonds, respectively. 
The beam FWHM is presented with gray hlled circle. 






53 



RA offset (arcsec) 


Fig. 19.— Integrated intensities of blue- and redshifted emission of CO 2-1. Blue contours 
represent intensities integrated over the blueshifted range from -20 to -7 km s“^ and red 
contours show intensities integrated over the redshifted range from -2 to 10 km s“^. The 
contours start from 4a of the peak intensity to 90 % and increase in steps of 10 %. VLA 6, 
VLA 7, and VLA 10 are denoted by crosses. The beam FWHM is presented with gray hlled 
ellipse. 



